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Abstract 



A theory is developed and applied to the study of opportunities and specific 
i-C I features of coherent control of four-wave mixing as well as of the accompanying 

processes in the continuous-wave regime, which involve transitions between bound 
^ . and free quantum states. Such opportunities become feasible through construc- 

ts I five and destructive interference of quantum pathways. Two coupling schemes of 

practical importance are investigated. In the first, a ladder energy level scheme, 
fully-resonant sum-frequency nonlinear-optical generation of short-wavelength ra- 
diation driven by several strong fields is investigated. The relaxation processes as 
^ ' well as absorption of the fundamental and generated radiations, which play an im- 

■ portant role, are taken into consideration. It is shown that the generation output 

can be considerably increased through the appropriate adjustment of several laser- 
induced continuum structures. In the second, a folded scheme, a possible control 
of two-photon dissociation (A-scheme) using auxiliary laser radiation applied to 
the adjacent bound-free transition (^-configuration) is investigated. Besides disso- 
ciation, the proposed method enables one to control population transfer between 
two upper discreet levels via the lower-energy dissociation continuum, while direct 
transition between these states is not allowed. The opportunities of manipulating 
these processes as well as of four-wave-mixing-based spectroscopy are explored both 
analytically and through a numerical simulation for Na2 dimers. 
PACS: 42.50.Hz, 42.50.Ct,42.65.Ky 42.65.Yj, 32.80.Fb., 32.80.Qk 
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1 Introduction 

The opportunities of laser control of optical properties and chemical reactivity of sub- 
stances associated with the continuum energy states of atoms and molecules, and based 
on nonlinear quantum coherence and interference processes, have attracted considerable 
attention since their theoretical prediction 0, 0] and experimental realization . These 
effects were shown to provide the feasibility of manipulating the spectral characteristics of 
the absorption, refraction, and nonlinear-optical conversion of weak radiations by an addi- 
tional strong field, which is in resonance with a transition between an excited vacant level 
and a state in the continuum. The strong-field effects in such coupling configuration man- 
ifest themselves by the appearance of spectral autoionizing-like laser-induced continuum 
structures (LICS). In the presence of a real resonance autoionizing level a strong field has a 
strong infiuence on the spectral characteristics of this resonance too 0, H, IE • The effects 
of the LICS on the absorption, refraction, nonlinear-optica l ge neration a nd p hotoioniza- 
tion may differ significantly (see, for example, Refs 0, B S Eo, IHl 13 0, E, [3] and 
the references therein). Most of the recent publications consider LICS in photoionization. 
The infiuence of real autoionizing levels and of the LICS on nonlinear susceptibility and 
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four-wave mixing (FWM) has been investigated for frequencies of the other fields far from 
one-photon resonances The interference of two LICS was shown to bring new 

quahtative effects in spectral properties of absorption, refraction and nonlinear optical 
generation of short- wavelength radiation [l6l |. 

Recently, the laser and chemical communities have focused on coherent quantum 
control of chemical reactions and other dissipatioii processes like photodissociation of 
molecules and photoionization of atoms JHS HBMB 

13, III H H 0, H H H H III H Ell 

In many cases, such control is based on implementation of the interference of different 
quantum channels, which give rise to laser-induced continuum structures (LICS) in the 
dissociation and ionization continua. Recent experiments aimed at control of photoion- 
ization of metastable helium atoms and some other processes making use of this technique 
are reported elsewhere 54, 5^ 56l |. 

This paper further develops theoretical approaches to these problems and reports novel 
opportunities for manipulating spectral characteristic of short-wavelength generation and 
photodissociation through the interplay of two LICS induced by strong continuous-wave 
laser fields coupling the dissociation continuum to two different bound states. Primary 
attention is paid to two important applications: (1) coherent control of sum-frequency 
short- wavelength generation in a ladder-type energy-level scheme (depicted in Fig. ^ and 
(2) coherent control of dissociation, population transfer and four-wave mixing feasible in 
a folded-type scheme (depicted in Fig. E)). 

As concerns the first problem, we generalize our investigation of the LICS to the case of 
several strong control fields in situations when one can expect the strongest enhancement 
of the generated power. We consider the combined influence of the interference of two 
laser-induced structures in the continuum and of the laser-induced transparency (for the 
one-photon-resonant initial radiation and for the generated radiation) on the processes of 
nonlinear-optical generation of short-wavelength radiation. We demonstrate for the first 
time that nonlinear interference resonances in the output power of the radiation generated 
in an absorbing medium may differ considerably from the corresponding resonances in a 
transparent medium because of the combined influence of the nonlinear resonances in the 
nonlinear polarization and in absorption of both the initial and generated radiations. 

A two-photon dissociation (A-scheme) controlled by an auxiliary laser field coupling 
continuum with another discrete state (^^-configuration) is investigated for the second 
problem. All the coupled radiations are strong enough to drive molecular transitions. The 
proposed technique enables the coherent laser control of dissociation related to the lower 
laying molecular 'dark' states that are not connected with the ground one by the allowed 
transition. On the other hand, the scheme under consideration enables one to transfer 
the population between two upper bound states which are not connected directly by the 
allowed transition. An analytical solution of the corresponding master density matrix 
equations of the problem is obtained, and a numerical analysis for relevant experimental 
conditions 0, |5y| is performed. The possibility of manipulating the dissociation spectra 
and populations of excited states at the expense of the interference of quantum pathways 
through a variety of continuum states is explored. The dependence of the effects on the 
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composed Fano parameters related to high excited levels, on the detuning between two 
Lies and on the intensities of the laser radiation is investigated. 

2 Four- wave mixing, absorption and refractive in- 
dices in ladder energy-level schemes controlled by 
interference of two LICS 

2.1 Equations for coupled electromagnetic waves and their so- 
lution 

Let us consider four plane-polarized electromagnetic waves travelling along the z axis of 
an isotropic medium, 

E\z,t) = Re{Ej{z) exp[i{ujjt - kjz)]}, (2.1) 

where kj is the complex wave-number corresponding to the frequency uj {j = 1, 2, 3, S). 
We assume that the fields Ei and Es are weak compared to the driving fields E2 and 
£"3, which do not vary along the medium. On the contrary, the fields Ei and Es can 
change considerably along the medium, because of both absorption and nonlinear-optical 
conversion. Then the spacial behavior of the waves Es and Ei is described by two coupled 
reduced wave equations, 

dEs{z)/dz = i2'iTk'sxP E2E3Ei{z) exp{iAkz), 

dEi{z)/dz = i2TTk[x? E;e;Es{z) exp{-iAkz). (2.2) 

Here kj = kj — ik'- = {2nujj/c)xj, k'- = aj/2, and XjyCij are the effective linear suscep- 

(3) (3) 

tibilities and absorption indices for the corresponding radiations, and Xi iXs 
nonlinear susceptibilities describing the four- wave mixing processes: us ^ uJi + uj^ + oj^, 
A/c = ks — ki — k2 — k^. The quantum conversion efficiency of the radiation Ei into Es 
varies along the medium as 

= {k[/k's)\Es{z)/E,{0)\'exp{-asz). (2.3) 

Let us first consider the case of low efficiency, for which the change in the Ei caused 
by the nonlinear-optical conversion can be ignored. Then the second equation in (j2.2p 
can be ignored as well and, with the account of the boundary condition Es{z = 0) = 0, 
one obtains the following for the generated radiation and quantum conversion efficiency: 

Esiz) = {27rkyAk)xfE^E2E,[exp{-iAkz) - 1], (2.4) 
r/q(2) = k'ski\\27rxfE2E3\'^/\Ak\'^)exp{-asz)\exp{-iAkz) - l|l (2.5) 

If the medium length is much shorter than the minimum absorption length Labs = 
min{Li = 2/ai,Ls = 2/as} and both of these are assumed much shorter than the 



4 



coherence length Lcoh = Afc' ^, then Eq. ()2.5p reduces to the approximate formula 



r]^ = ks'h'\2nxfE2EsrL 



|2r2 



(2.6) 



where Lg represents either the length of the medium (in the case of weak absorption) or 
the optimal length of the order of m.m{Labs, Lcoh}- 

For a medium with substantial absorption dispersion (ai 7^ as), but Ak' = and 



(3) 

Xi 



(3)* 



Xs' , the solution of equations ()2.2|) . for the more general case of considerable 
conversion, takes the form 



Vciiz) = 4— exp 



[ai + as)- 



sinh 




, (2,7) 



Here figo = k[k's I 2tixs E2E^ \' is the conversion efficiency per unit of the medium length 
under constant fundamental radiations, and h = AfjqQ — [ai —05)^/4 defines the difference 
between the rates of absorption and nonlinear-optical conversion of the radiations. If 
6 < 0, the conversion rate of hwi photons into ^5 photons is less than their absorption 
rate; if 6 = 0, the photon absorption and conversion rates are equal; and if 6 > 0, the 
nonlinear-optical conversion rate exceeds the absorption rate. In the latter case, one can 
expect an oscillatory dependence of the transfer of the weak radiations from one to the 
another and back along the medium. 

The susceptibilities X^s^ ^ Xi \ Xs ^-nd Xi can be derived from the medium polarizations 
at the corresponding frequencies: 



Piujj)=XjE,, P^L(^5) 



xfE^E^E^, 



p^L(^^) = ^^^^Ese;ei 



(3) 



(2i 



These components can be calculated conveniently with the aid of a density matrix, pjj. 



as 



V{u}j) = Npij{ujj)dji + c.c. 



(2.9) 



where is the atomic number density in the medium, and dji is a matrix element of 
the projection of the transition electric-dipole moment along the direction of the electric 
vector of the corresponding field. The problem of finding and optimizing the quantum 
efficiency of the conversion process thus reduces to finding the off-diagonal elements of 
the density matrix. 



2.2 Density matrix master equations and their solutions 

First, we shall consider the transition scheme depicted in Fig. ^ A strong field E2 at 
frequency UJ2 is close to resonance with the transition between levels m and n, while the 
strong fields E and E3 at frequencies uj and us are in resonance with the transitions 
between levels / and n and some states in the continuum. Radiation at the frequency us 
can be either a probe or generated by four-wave mixing. We shall investigate the infiuence 
of the strong fields on the spectral characteristics of the absorption of the radiations Ei and 
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Es at the frequencies uji and us (both as independent probe radiations or, alternatively, 
as the frequencies linked through generation, us = UJ1+UJ2 + 1^3)- The field Ei is assumed 
close to resonance with the transition from the ground state to the level m and Es to 
the transition from the ground state to the continuum. These radiations are weak, so 
that a change in the level populations due to all resonant couplings can be ignored. We 
neglect the degeneration of all states, including those in the continuum. In addition to 
the contribution of the off-resonant continuum states, we account for the contribution of 
discrete off-resonant levels which are combined to form the level k in Fig. ^ We suppose 



that the detunings | ui-ujgm 



UJl + UJ2 



UJ. 



gn 



and 



UJ3 



UJ. 



nf 



are considerably 



less than all the other detunings. The equations for the density matrix, considered in the 
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Figure 1: LICS-based coherent control in ladder schemes. 

interaction representation to within the first order of perturbation theory in weak fields, 
can be written as follows: 

gm/ dt -\- ^ gmPgm ^(^Pgg^gm ~l~ Pgn^nm) y 
dpgn/ dt + TgnPgn = i J PgeVends + iPgrnVmn + 1 PgkVkn, 

dp,f/dt + r,fP„ = .jp,,%sde + ij:p,M„ (2,10) 

dpge/dt = i{PggVge+PgnVne+PgfVfe), 
dpgk/dt + TgkPgk = KpgnVnk + PgfVfk)- 

Here, the index e denotes the continuum states; Vmn = GmnGxp[i{uj2 — uJnm)t],Vgs = 

Ggeexp[i{uJs - LJgs)t],Vns = Gneexp[i(u;3 - LJne)t],Vgm = Ggmexp[i(cJi - LJgm)t],Vfe = 

Gfeexp[i{Lj-Ljf^)t],Vkn = G kn Gxp[i{uj - uJkn)t], Vkf = Gkf exp[i{LJ3 - LJkf)t] are the matrix 
elements of the Hermitian interaction Hamiltoman V, considered in the electric-dipole 
approximation and in the interaction representation (in units of h); Vij = V*f,Gmn = 
—E2dmn/'2KGg£ = —E4^dg!r/2h, etc.; and Tij is the homogeneous half-width of the i — j 
transition. In the approximation of the weak fields Ei and Es, we obtain pgg ^ l,pm ~ 

Pn ^ pf « 1. 
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Under steady-state conditions, each off-diagonal element of the density matrix is a 
sum of two components, which may oscillate at different frequencies: 

Pge = rge exp[i(u;5 - ujg^)t] + Rg^ exp[i(u;i + lj2 + uj^ - ujg,)t\, 

Pgn = Tgn exp[i(t^5 - UJ^ - UJgn)t\ + Rgn exp[i(u;i + UJ2 - UJgn)t], 

Pgm = rgmexp[i{uJs -UJS-UJ2- UJgm)t] + exp [i(u;i - UJg,n)t\, (2.11) 

Pgf = ^gf'^^PK^s -Ld - ujgf)t] + i?g/exp[i(u;i + uj2 + UJ3 - lj - ujgf)t], 
Pgk = Tgk exp[i(u;5 -uj^-uj - iOgk)t] + Rgk exp[i(u;i +uj2-uj - idgk)t\. 

By substituting ()2.1H) into ()2.10|) . one can see that the set of differential equations under 
consideration reduces to two independent sets of algebraic equations, where each refers 
to the processes determined by only one weak field: 

'^RgmDgui G gui Rgn.Gnmi Dgui ^ gm ~\~ i(^^l ^grri)i 

■^RgnDg^ = -j RgsG.r4e - RgmG„.r. " Z^,„ = T g^ + i^U, + U2 - Ug.), 

'^RgeDge = —RgnGne " RgfGfe, Dg^ = i{uJi + UJ2 + ^3 — ^ge), (2.12) 

iRgkDgk = —{RgnGnk — RgfGfk), Dgk = Tgfc + ^{^^1 + UJ2 — UJ — Ulgk), 

irgePge = -Gge - TgnGne - fgfGfe, Pge = i{^S - ^ge), 

^^gmPgm — ~^gnGnm, Pgm — ^ gm + K^S — UJ3 — UJ2 — ^gm)-, (2-13) 

ir,/P»/ = - / r„G,,de - ,,,G,„ p„ = T,, + - . - 

irgkPgk = -TgnGnk " ^gfGfk, Pgk = gk + i{^s - uj3-^ - ^gk)- 

Here and later the repeated index k indicates summation over all discrete off-resonant 
levels combined to form the level k. 

Equations ()2.12|1 describe the absorption of Ei and generation at the frequency ujs, 
whereas ()2.13|) presents the absorption of Es and parametric conversion of Es back into 
El. One can solve (12) by substituting the third equation into the second and fourth. 
Then, in the calculation of the integrals, one can employ the (^-function, 

[-l{uJs - UJeg)]~^ = 7r5{u;s - UJ^g) + i'P(c^5 - ^eg)~^, (2.14) 

where 5{^) is the delta function, and V is the principal value obtained by integration. 
This leads to the following equations: 



gei 



Rge = i[Gne - Gfe{lnfhff)f3f{l - IQn f) / X f]Rg J D 

■^Ggm ~l~ RgnGnm ^ GgmGuin-Xf 
^ Xm^gm ' ^ ^gf^gnXmi^ + gnn){XnXf — K + AmXf)' 
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where 



(^f = 9ff/i^ + 9ff), Pn = 9nn/{l + 9nn), 9nn = Inn/^ 
K = (3f(3n{'\- — ^QnfY, Am = gmn/ Xmil + gnn), Qmn 
Qij = ^ijhij, lij = T^^GieGej \e=nws ikG kj / P gk) ■, 



9"' 



I a 



9ff = Iff/^gf^ 

P/r r • 



6.. 



hP 



de- 



lm[GikGkj/P: 



'gk\ 



= 1 + iXi, Xn = (fil + Q2 — Snn)/(Xgn + Inn), Xm = Vti/Tgm, 
Xf = (l^i + ^2 - l^L - '^//)/(rg/ + 7//) = {U1+UJ2 + UJ3-U- Ugf - Sff)/{Tgf + 'Jff), 
= {UJ + UJgf) - {uJs + UJgn), = UJl - UJgm, Vt2=UJ2-UJmn- (2-16) 



Here fi^ is the spacing between the quasi- levels induced by the radiations E and in 
the continuum. The Fano parameters qij [s^ are assumed real and indicate the ratio 
of the light-induced shifts and the broadening of the corresponding resonances by the 
control fields. In the adopted approximation, these parameters are independent of the 
field intensities and are governed solely by the properties of the investigated atom. 

Following the same procedure as above and bearing in mind the condition 0;^ = 
uji + UJ2 + uJ^i one finds from the set of equations ()2.13|) that 



r, 
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l{Gge - Gfe{.lgfhff)Pf{l " %/)/^/ + rgn[Gne " G fe{.lnf /l f f) 13 f{l - IQuf) / X f]} / Pg,, 
_ (1 - ig/n)(l - '^qgfhgflfn/iXgf + 7//) - (1 - igg„)X/7g„ 



r G 

. I qn^nm 
I gm — ^ ^ -p 1 ' gn 



Y V 

-^m^ gm 



gnn)'- gn 



(2.17) 



The calculation of Rgm from ()2.15|) and r^^ from ()2.17|) . and application of formulas ()2.8|) 
and ()2.9|) after integration over the continuum states, gives the following expressions for 
the absorption and refractive indices at the frequencies uji and a;^, respectively: 



aio 

"so 
Fs 



ReFi 
1 



n{iJi) 



2(^1 max - 1) 

A-rnXf 



ImFi, 



Xr, 



ReFs 



1 - 



XnXf + AraXf 

n{u)s) - 1 



2 (ris max 
XfXn{A,r, 



A 



f) 



K 
ImFs, 

- U 



(2.18) 



AmAfXf 



Af-A, 



1-Af-An- 



XfXn — K + AmXf 

K{A^ + Af)-U-AmA^Xf 



XfXn — K + AmXf 

K{An + Af)-U ' 
XfXn - K . 



(2.19) 



Here ctio is the absorption index with the fields intensities and detunings turned to zero, 
and a^o is the similar quantity for a transition to the continuum; riimax and n^max are 
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the maximum values of the corresponding refractive indices under control fields turned 
off; and 

Af = f3f{l - iqgff/Xf, An = - iqgnf/Xn, U = 2/5//5„(l - iqgf){l - ig/„)(l - ig„g). 

The functions A^ = — iqgnY / and X^ = X^ + A„i account for the perturbation 
of a two-photon resonance with the level n by the strong fields. The expressions for the 
refractive index are obtained on the assumption that this index is close to unity in the 
absence of the fields. 

The calculation of Rge and r^m from the expressions (j2.15|) and (j2.17p . which is car- 
ried out making use of formulas (j2.8|) and ()2.9|) and integration over the states in the 
continuum, gives expressions for the FWM nonlinear susceptibility at us = uji + uj2 + co's 
as 

X^'K^s) ^ Xf-Pj{l-iq^f){l-\qjg)/{l-\qng) 

^g) Xm{l+gnn){XnXf-K + A^Xf) ' > 

and for the nonlinear susceptibility determining conversion of the Es radiation back into 
the radiation of frequency uoi as 

X^^KuJi) _ Xf- (3f{l - iqnf){l - ig/,)/(l - iqng) , . 

^(3) Xn,{l + gnn){XnXf-K + AmXf) ' 

In the expressions ()2.20p and ()2.21|) . the quantities X50 and Xio^ are fully- resonant non- 
linear susceptibilities for non-perturbative weak fields. 

2.3 Laser- induced structures in absorption and refraction 
spectra at discrete and continuous transitions 

In this subsection, we analyze and compare the effects of strong fields on discrete and 
continuous spectra with the aid of the derived formulas. An important difference is 
seen as compared with similar effects if solely discrete transitions are involved [sil . l6ol | . 
The results given below demonstrate considerable perturbations of discrete spectra by 
the radiations coupled to the continuum. The first term in ()2.18p represents the field- 
unperturbed absorption coefficient for the gn transition, and the second term refers to 
the cumulative effects of the strong fields. The coefficient Am. represents splitting of 
a resonance into two components by the strong field E2 jsil l6ol | and modified by the 
strong field E^. The function K describes modification of these effects by the strong 
field E. It is proportional to the product of the intensities of the fields E and E^. The 
effects in question disappear when any of these fields is turned off. Since the field E2 
is in resonance with a discrete transition and the fields E^ and E are coupled to the 
continuum, the spectral properties of the corresponding contributions are different. If 
-E'3 = (7„ = 0,/3„ = 0,i^ = Q.gnn = 0), the equation for absorption in ()2.18|) converges 
to the standard one for three-level nonlinear spectroscopy js^ 60], 



a(ci;i)/aoi = Re 



^ + gn 



(2.22) 
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Figure 2: Absorption index at Ui (reduced by its resonant value in the absence of all strong 
fields) vs one-photon detunings Qi/Tgm- Here, qjf = 0.9, g„„ = 0.5, Tgm/Tgj = 100, and 
^gmf^gn = 10. (a,b): fig = 0, Qfn = 1.5, gmn = 7, Iff/Tgf = 2. (a): 7„„ = 0(1), 

Innf^gn = 5, ^Lf^gm. = 0.8(2). (b): -fnuf^gn = 5, fi^ = 0. (c,d): n2/Tmn = 0.3, 
9mn = 70, Iff/Tgf = 10, fii/rgm = -1.1, 'Jnn = 0(1), 7„„/rg„ = 50(2). (c): Qfn = 15. 

(d): Qfn = 1.5. 

The factor in the parentheses in the above formula has two singularities with respect 
to Xm, which indicates splitting of the resonance into two maxima. The positions of these 
maxima and their relative amplitudes may vary depending on the parameters of the fields 
and transitions. Resonance splitting is stipulated by the appearance of coherence at the 
transition ng and by the consequent appearance of additional quantum transitions in 
which photons of frequency hui may participate. The phase and relaxation properties of 
the corresponding term in nonlinear polarization at the frequency Ui are represented by 
the dispersion function X„. The additional strong fields and E perturb the quantum 
system, which, as pointed out earlier, causes additional modification of the laser-induced 
quasi-levels and resonances. The changes in the refractive index can he explained similarly. 

The spectral characteristics of absorption at the frequency us are governed by the 
interference between three quantum pathways: directly to the continuum, to level /, and 
also to the superposition of levels n and m. The influence of the first two processes was 
investigated earlier in Refs. [l|, 0, H, Ig^] and is accounted for by the function Ai in the 
equations (j2.19|) . The strong fields E^ and E2 lead to additional changes in the spectra. 
The additional independent structure, described in these expressions by the function An, 
is supplemented by the interference structures proportional to the functions K and U, 
which disappear when any of the fields E3 or E is turned off or when the spacing between 



10 



the quasi- levels is increased. 

Figure |21 shows the dependence of the absorption index at uJi on the scaled detunings 
from the bare-state one-photon resonance VLi/Vg^- The dashed plot at Fig. |2fb) shows the 
lineshape of the dispersive index at the frequency uji for the same parameters. The plot 
indicates very strong dispersion, induced by the dressing fields coupled to both discrete 
and continuum states. Figures Efc) and Efd) illustrate the sensitivity of the spectra, 
induced jointly by two dressing fields and on the Fano parameter One can 
create transparency in certain frequency intervals of the discrete transitions or, on the 
contrary, eliminate effects of the dressing fields with the aid of destructive interference by 
varying the intensities and detunings of the driving fields as well as the coupled levels. 
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Figure 3: Absorption index at ujs (reduced by its value in the absence of all strong 
fields) vs detunings f^/F/g = {ojs — uj — ujfg)/V fg (a-c). The inserts are the interference 
contributions to the corresponding curves vs detuning (for (b,c) the detuning interval is 
the same as for the main curves). Here Vg^/'^gf = 100, Tg^/Vgn = 10. (a,b): qjf = 0.9, 

Qnn = 0.9, ^Lf^gf = "HO, and ^2/'^gf = 30. (a): Qgf = -0.5, Qgn = -0.95, Qfn = 15. 

(b): qgf = -0.95, qgn = -0.5, g/„ = 15(1,2,3), g/„ = 150(4). (1) ^2 = Eg = 0, 
jff/Tgf = 10. (2-4) gmn = 70, 7nn/r,. = 50. (2) 7// = 0. (3,4) ^ff/Tgf = 10. (c): 
Qgf = 0.95, qgn = 0.01, qff = 0.01, g„„ = -5, g/„ = 1.5, ^^2/'^gf = 0, Iff/Tgf = 10, 
F,^ = -1530. (2,3) 7nn/r,„ = 5. (2) Q^rgf = -110. (3) fi^/F,; = -405. (d-f): 
Scaled dispersive profile of LICS at us vs fl/Tfg. (All parameters are the same as in plots 
(a-c), consequently.) In the absence of the driving fields absorption level corresponds to 
1, dispersive - to zero. 
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Figures Efa) and ^h) show the dependence of the absorption index at the frequency 
ujs (reduced by its value in the absence of the dressing fields) on the scaled detuning 
O/F/g = {us — oj — ujfg)/Vgf. The plots indicate the sensitivity of the probe field absorp- 
tion on the Fano parameters (in the considered cases on qgf and qgn) and demonstrate 
new possibilities for manipulation of LICS with extra dressing fields coupling adjacent 
transitions. Interference contributions, which disappear in the absence of either the E or 
£^3 field, are shown in the inputs to the figure. Figure Efc) demonstrates the sensitivity of 
the absorption spectrum on the strength of the driving field -Es as well as on the spacing of 
two quasi-levels, induced in the continuum (on two-photon detuning VLl = ^ — ^3 — ^fn)■ 
^he figures show possible manipulation of LICS in the absorption index for a probe ra- 
diation (including formation of transparency windows). This is brought about by the 
interference of two LICS (quantum pathways via discreet and continuum states), induced 
by the fields at and uj and modified by the strong field at uj2- 

The refractive index is proportional to the derivative of the structures over frequency. 
The plots infd,e,f) display the potential of manipulating the magnitude and lineshape of 
the factor {n{ujs) — l)/(?^5max ~ 1) with the change of intensities and detunings of the 
dressing fields, as well as strong dependence on the Fano parameters, which may find 
applications in short-wavelength optics. 

The nonlinear interference effects involving quantum transitions may influence dif- 
ferently the linear and nonlinear susceptibilities. Therefore, under certain conditions a 
reduction in the absorption of the initial and generated radiations can be simultaneously 
followed by an increase in the nonlinear polarization caused by the resonant effects and 
by the constructive interference in the strong fields. It is thus possible to increase consid- 
erably the power of the generated short-wavelength radiation. 

2.4 Resonance sum- frequency generation in strongly-absorbing 
media enhanced by quantum interference 

In weak fields the nonlinear susceptibility increases strongly upon approaching to discrete 
resonances, but this is accompanied by enhancement of the absorption of the initial ra- 
diations. Depending on the detunings from the resonances, on the ratio of the oscillator 
strengths of the transitions, and on the radiation intensities, either the absorption of the 
initial radiation Ei or of the generated radiation may predominate. A numerical analysis 
of the influence of these factors on changes in the frequency dependence of the power of 
the generated short-wavelength radiation during the course of propagation in an optically 
dense medium can be made if the quantum efficiency of conversion given by expression 
fl2.7|) is rewritten in the form: 

rj^(z) = ^exp[-{a,+Cas)zo]!^smh' .J {\h\ -h)C {\h\ +h)C |, 

(2.23) 

where f/qo = ^qo/("io"5o), & = &/("ioaso) = 4r/qo - («i - Ca5)V4C, ai = ai/aio; as = 
Ois/oiso, C = a^o/ttio, andzo = zaio/2. 
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The expression for rj^Q, describing the quantum efficiency of conversion over a distance 
I/a/choQ^so; considered within the approximation of ignoring absorption of given fields, 
becomes 

VciO = VqolX^^^ \^9mngnn, (2.24) 

where rj^Q is the quantum efficiency of conversion based on the resonant unperturbed 
nonhnearity over a distance l/^aio«5o in fields corresponding to gmn = 9nn = 1; X^'^^ = 
X^^Vxo- We shall henceforth use the following approximate expressions: 



aio = 4:7iLJi\dgm\^/chTgm, aso = 4:7i'^{us/c)\d 



gel 



\Xo\ = {TT/2h) {l + qgj\dgmdmndnedeg\ C^gm^gn) ■ (2.25) 

In this approximation the factor ?7qQ is determined completely by the Fano parameter g^^: 

r/Jo = k[k's\27iXo\^\/d^r.dne\M^Qfi^^grnTlJ-'a^oaso = (1 + ^ J- (2.26) 

As pointed our earlier, the laser-induced absorption and transparency resonances of 
the radiation Ei can be interpreted as splitting of the level m into quasi-levels by the 
strong field E2 (in combination with the fields E3 and E). The laser-induced resonances 
of the generated radiation are determined by two quasi-levels in the continuum and appear 
near the frequencies Ung + uj^ and Ufg + These quasi-levels are separated by an energy 
HQl- The detuning of the generated radiation frequency from the ffist resonance amounts 
to ^5 = f2i +f22, and that from the second resonance isf2 = r2s' — ^2^. The corresponding 
quantum channels may interfere differently when the detunings Qi, Q2 and are altered, 
which is manifested in the spectral characteristics of the absorption and nonlinear-optical 
generation processes. The relative role of these channels is governed by the intensities of 
the radiations, by their detunings from the resonances, and by the ratio of the oscillator 
strengths for the g — m and g — e transitions. We shall now illustrate these relationships 
by considering several numerical models. 

Figure m illustrates the case when the integral oscillator strength for a transition to an 
energy interval E (of the order of the level widths) in the continuum is considerably less 
than for the gm transition {C = 10~^). The frequencies Ui and 002 are detuned far from 
their one-photon transitions, but the frequency sum is close to a perturbed two-photon 
resonance. It is evident from Fig. |3]that changes in Ql (because of variation of the fre- 
quencies u or CJ3) reduce, for the selected Fano parameters and radiations, the absorption 
coefficient ai by approximately threefold, increase as approximately by a factor of 3.8 
in one detuning interval, and reduce considerably this coefficient in the other interval, 
whereas the square of the modulus of the nonlinear polarization (proportional to r/qp) 
increases by a factor of 1.9 (these changes are relative to the values of the corresponding 
parameters in the far wings where the effects of the strong field E do not appear) [Fig. 
131(a)]. The absorption coefficient for a transition to the continuum remains on the whole 
much less than ai (a^/ai = Cas/'^i ~ 10~^). In a certain range of fi^, the sign of b 
becomes positive, and the nonlinear-optical conversion rate begins to exceed the rate of 
absorption of the radiation. It is in this interval that there is a sharp maximum of the 
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Figure 4: Dependence of the absorption indices ai / aio and as / aio, of the value rj^Q (that is 
proportional to the squared modulus of the amplitude of the nonlinear polarization at the 
entrance to the medium) , and of the conversion rate h on the detuning VL l (a) . Dependence 
of the quantum conversion efficiency ?7q on VLl for zaiQ = 8.5 ■ 10^ (1), zaio = 1 ■ (2), 
and zaio = 2 ■ 10^ (3) (b). Dependence of the quantum conversion efficiency ?7q on the 
optical thickness and on Ql (c). Here, C = 10~^, gff = 150, gnn = 200, gmn = 9000, 
ni/Tgf = 5000, n2/Tgf = -5100, Qfg = 0.95, qgn = -2, qff = 0.01, g„„ = -5, g/„ = 0, 

TgjTgJ = 100, TgjTg^ = 10 ( C ) . 

quantum efficiency of conversion (the dependence of this conversion efficiency on the op- 
tical thickness and on Ql is demonstrated in Fig. Il](c), where the maximum ?7qmax = 0.29 
is reached for zaio = 4000). The absorption of the radiation Ei alters considerably the 
dependence of the power of the generated radiation onQi along the medium [Fig. Il](b)]. 

Figure El is computed for the case where the detuning from a one-photon resonance 
is still large, but the oscillator strengths for a discrete transition and a transition to the 
continuum differ less (C = 3 - 10~^). Then, under the action of the field E, variation of Ql 
can reduce the absorption coefficient ai by a factor of just 1.5, whereas the absorption co- 
efficient as increases approximately threefold in one interval, but falls considerably in the 
other interval; practically throughout the whole range of the detuning Ql, the dominant 
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Figure 5: Absorption indices ai/aio and ag/aio, the value r/qp and conversion rate b 
plotted as a function of Ql (a). Dependence of the quantum conversion efficiency rj^ 
on Ql along the medium computed for zaio = 4.5 ■ 10^ (1), zaio = 5 ■ 10^ (2), and 
zaio = 5.5 ■ 10^ (3) (b). Dependence of r/q on the optical thickness and on Ql (c). Here, 
C = 3 ■ 10"^, gnn = 500, gmn = 8000 (a-c). The other parameters are the same as in Fig. 

H 

effect is the absorption transition into the continuum {as/di ~ 70), and TJ^q increases 
by a factor of 1.9 [Fig. Efa)]. The quantity b is always positive and has a characteristic 
maximum in a certain interval of f^i,. For a given detuning, the rate of conversion begins 
to exceed the rate of absorption so much that an oscillatory regime appears along the 
medium. Propagation in an optically-dense medium makes the spectrum of the output 
power significantly dependent on the optical length of the medium, i.e., on z or on the 
concentration of atoms. The absolute amplitude and the position of the maximum 
change, and new resonances appear [Fig. Efb)]. At the first maximum (corresponding to 
zaio ~ 125) the quantum efficiency of conversion can reach 0.9 [Fig. Efc)], which is three 
times greater than in the preceding case. 

In view of the nonresonant nature of the interaction, the attainment of such high val- 
ues of the quantum efficiency of conversion requires high intensities of the E2 and E^, 
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Figure 6: Dependence of the absorption indices ai/aio and a^/aio, of the value tj^q, and 
of the conversion rate b on Ql (a). Dependence of the quantum conversion efficiency r/q 
on the optical thickness for fii/r^/ = (1), fii/Fg/ = -250 (2), flL/Tgf = -400 (3) (b). 
Dependence of the quantum conversion efficiency rjq on the optical thickness and on 
(c). Here, gfj = 100, gnn = 5, gmn = 7, ^i/Tgj = 0, ^2^^^/ = -250 (a-c). The other 
parameters are the same as in Fig. El 

radiations, and also long lengths (or high densities of atoms) of the medium. The use 
of resonant processes makes it possible to reduce the required intensities and to reach 
considerable values of the quantum efficiency by optimization of the bleaching and inter- 
ference effects. Figure El illustrates the case where the ratio of oscillator strengths is the 
same as in Fig. but there is no detuning from a one-photon resonance. At intensities of 
the E2 radiation three orders of magnitude less than in the preceding case, it is possible 
to reduce the absorption coefficient of the Ei radiation approximately by a factor of 10 
compared with the value of this coefficient in the absence of strong fields (the maximum 
effect of the field E is a reduction in this coefficient by a factor of 1.5) [Fig. Eta)]. In 
the region of reduction in ai, the value of as increases approximately threefold and r/qo 
by a factor of 4.7. The absorption coefficients ai and as are then comparable, and the 
nonlinear-optical conversion rate exceeds considerably, in a wide range of Ql, the rates of 
absorption of the radiations, reaching a sharp maximum at Q^/Ygf = —250. This gives 
rise to an oscillatory regime of generation along the medium in this detuning range [Figs. 
Efb,c)], i.e., it leads to the feasibility of total (apart from that lost by absorption) con- 
version of the radiations Ei and Es (and vice versa) for certain products of the length of 
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the medium and the concentration of atoms. For the selected parameters, the conversion 
efficiency at the first maximum (where zoio ~ 5) reaches 0.54. This is less than in the 
preceding case because of the much greater absolute values of the absorption coefficient, 
but this efficiency is reached at much lower intensities and for much smaller lengths of 
the nonlinear medium and much lower concentration of atoms in this medium. 

2.5 Absorption and dispersive spectra at Doppler-broadened 
transitions 





Figure 7: Absorption (solid) and refractive (dash) indices at ui (reduced by its value in 
the absence of all strong fields) in a Doppler-broadened medium vs detunings Qi/Aum- 
Here, the Doppler HWHM Auid = IQ.GbTgm, the wavevector orientations are k ff ki, 
k2 TT ks II ki, and k2/ki = 0.9, k^/ki = 0.5, k/ki = 0.6 (kj is wave-vector corresponding 
to the frequency uji), Tg„jTgf = 100, Tgm/Tgn = 10; \Gr,^n\'^ / {AuiDy = 1, g„„ = 0.5, 
Qff = 0.9, n2/AuiD = 9. (a,b): 7„„/Au;iz) = 0.2, 7///Au;id = 0.1, g/„ = 0.5. (a): 
Ql = 0. (b-d): Ql/Auid = -0.8. (c): -fnn/AuiD = 0.8, 7///Ac^id = 0.3, g/„ = -1.5. 

(d): 'Jnn/AuJiD = 0.2, -fff/AiOiD = 0.8, qjn = 1.5. 

A Lorentzian line profile with near natural linewidth can be observed only by making 
use special techniques, like atomic jets. Analysis shows that the contributions of atoms at 
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Figure 8: Absorption (solid) and refractive (dash) indices at ojs (reduced by its value in 
the absence of all strong fields) at Doppler-broadened transitions vs detunings VL/ Auj^^d = 
{los — (jJ — ujfg)/ AusD- Here, the Doppler HWHM Aljsd = 5 ■ lO^Fg/; the wavevector 
orientations are k || ks |t k2 tt ^4, and k/k4 = 0.8, A;3/A;4 = 0.3, k2/k4 = 0.37, 

r^m/r^/ = 100, TgjTgr, = 10; |G„„|V(AcUsd)' = 1, Inn/ACJSD = 0.4, ^ff/AcUsD = 0.8, 
Qgf = 0.95, Qgn = 0.01, Qff = 0.01, Qnn = "5. (a): Qfn = 1.5, fii/At^sz) = 1.5, 

02/Au;5D = 2.2. (b-d): 1^2 = 0. (b): = 15, Ql/Ausd = 1.5. (c): = -1.5, 
fii/AcusD = 15. (d): qfn = -1.5, Ql/Ausd = -0.5. 



different velocities in a Doppler-broadened medium may completely change the features 
described above. Figures [7| and |H1 demonstrate that proper adjustment of the orientation 
of the wavevectors, along with the intensities and detunings of the coupled waves, provides 
additional means to manipulate the lineshape. Moreover, enhanced subDoppler structures 
can be formed by the compensation for Doppler shifts by power shifts. For details of the 
physics, see 0, Ig^I and references therein. The plots dshow that the interference of 
contributions from the atoms at different velocities brings an important distinction in 
appearance of quantum interference processes at coupled discrete and continuous states. 
The figures show that the appearance of quantum interference at Doppler-broadened 
transitions also may be both constructive and distractive, depending on the detunings of 
UJ—UJ3 from Unf, on the Fano parameters, on the detunings from the two-photon resonance 
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gn, on the ratios of the wavenumbers, and on orientations of the wavevectors. 



3 Four- wave mixing, dissociation and population trans- 
fer controlled by LICS in folded energy-level schemes 

Here, we will consider the folded scheme, which is characteristic for molecules with low- 
lying predissociation states. Such a coupling configuration allows coherent laser control 
of dissociation through the 'dark' states that are not connected with the ground one by 
the allowed transition. On the other hand, the scheme under consideration enables one to 
transfer the population between two upper bound states which are not connected directly 
too. Thus, we will explore the possibilities of manipulating dissociation and populations 
of excited states through implementation of the interference of quantum pathways via a 
variety of low-lying continuum states. 




m 



Figure 9: LICS-based coherent control in the folded schemes. 

The proposed coupling scheme is illustrated in the energy level diagram depicted in 
Fig. ini The radiation at frequency ui couples the bound-bound transition m — n, and 
radiations at the frequencies UJ2 and couple the bound states n and / with the states of 
dissociation continuum e, as shown in the picture. Density matrix equations, which are 
employed for investigating population transfer, account for both relaxation and incoherent 
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excitation of the discrete states: 



Pne = i {pnnVne + PnfVfe) , 
Pfe = « (Pf/yfe + PfnVne) , 
Pme ^ iPmnVne ~l~ Pmf^fe) i 

Pmf + r^/Pm/ = i {^j dep„,eVef - V„,nPnf^ , (3.27) 
Pnf + ^nfPnf = depneVef - j deVnePeJ " KmPm/^ , 



Pnn ~l~ ^npnn 2Re ^ dEVnePen ~l~ i^nmpmn j ~\~ Qny 

Pff + T/P// = -2Re 1 rfe\//.p./^ + Q/, (3.28) 

Pmm ~l~ rmPmm 2Re (^iVmnPnm) ~l~ Qm ~l~ ^nmPnn- 

Here, Qj is the rate of incoherent excitation to level z, Tipa is the rate of relaxation, Wnm 
is the probability of excitation of level m due to the relaxation transitions from level n, 
and Vjk and Vj^ are the corresponding matrix elements of the quasi-resonant electrodipole 
interaction of the driving field and molecules scaled to Planck's constant h. The system of 
equations ()3.28p corresponds to the case of an open energy level configuration. This implies 
that all the levels coupled with the driving fields, including level m, can be incoherently 
excited from the reservoir of much greater populated lower levels, so that the parameters 
Qj are constant. 

For the closed schemes, i.e., the schemes where level m is a ground one, and therefore 
incoherent excitation of the upper levels depends on population transfer stimulated by 
the applied fields, the equations for the populations take the form 



Pff + TfPff = -2Re J deVfePef j + Wfpm 

Pnn + r„,p„„, = -2Re j deVnePsn + ^KmPmn ) + W 

npmm) 

(3.29) 

Pmm = I - pnn - Pff - I Wdt. 



Here, Wi describes the probability of excitation to level i from the ground state, and W 
is rate of dissociation. 
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In the equations ()3.27|) . the terms hke PekVkj (j, k = n,m) are discarded. This is vahd 
because coherence at discrete-discrete two-photon transitions is usually much stronger 
than that at bound-free two-photon transitions. The driving fields couple the interval of 
continuum states in the vicinity of the resonant energy Eq that is roughly equal to the 
width of the coupled power-broadened discrete levels. However, the oscillation strength 
at a bound-free transition is distributed over a much wider energy interval. Therefore, the 
fraction of oscillator strength attributed to the coupled energy interval is usually relatively 
small. The requirement derived with the aid of the Laplace transformation indicates that 
the indicated terms in the above equations can be neglected until 

df/dee=e, » h\V,k\d'f/del^,^, (3.30) 

where f{e) is the energy density of the transition oscillator strength. Such a requirement 
is fulfilled practically for all realistic atomic and molecular continua in the energy range 
that is well above the ionization or dissociation threshold. 



3.1 Quasistationary solution of density matrix equations 

In this section we shall assume that all radiations are continuous waves or rectangular 
pulses with duration r that is much greater than all relaxation times in the quantum 
system under investigation. We will consider also the case of small loss of molecules due 
to dissociation for the period of measurement or the pulse duration, so that the quasi- 
stationary regime is established. The corresponding requirement is 



T 

I 



Wdt = W{t) « 1, (3.31) 

where the dissociation yield Wdt is described by the equation 

W = p, = -2Re (^i j deV,nPne + i j deVefPf^ . (3.32) 

Quasi-stationary solutions of the density matrix can be found in the form 

Pa =ri [i = m, n, /), pij = rtj exp{iQijt), pj^ = rj^ exp{iQjet), 

Vij = Gijexp{iQijt), Vjs = Gjeexp{iQjet). (3.33) 

Here, Gij = Ekdij/2h, Gje = Ekdji;/2h, flmn = tOi — tOmn, ^nf = Us — U2 — Unf, ^mf = 
LUi - UJ2 + UJs - UJrnf, ^ne = ^^2 - t^ne, ^fe = ^^3 - ^^/e, and fi^e = UJi - UJ2 - ^me- The 

amplitudes of the resonant fields Ek (see Fig. IH)), and consequently all amphtudes r^, r^j, 
and r,£, are assumed to be independent of time. Then, by introducing the expressions 
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()3.33|) to the equations ()3.27|) . we can reduce them to a system of algebraic ones: 

irnePne = -rnGne " '^n/C/e, Pne = ^{.^2 " ^ne)] 
iffePfe = -rjGf, - TfnGne, Pfe = «(^^3 " ^fe)] 
irmePme = -TmnGne " ^mfG fe, Pme = ^(^^1 " " i^me)] 

.^.,^-j.r...-r^.,.0.^.^ ..^r,,..,.-..,. ,3.3, 

irmfPmf = - j dermeGef + Gmnrnf, Pmf = T-m/ + i{u;i - U2 + UJ-^ - Umf)] 

ir.sP., = - / der„G„ + / d.G„r„ + G„„r„.„ = r,„ ± ,(.3 - - 

Further by introducing the solutions of the first, second and third equations to the integrals 
and with the aid of the (^-function fj2.14p . we obtain 

"^mn \^nn ^mn. "I" '^^mn + 7nn.)] ~ '^m) Gmn '^mf {^fn ~^ ^'^ fn) i 

rmf [6fj - + l{V.^f + 7™ )] = GranTnf ' Tran (C/ + ^Inf) , (3-35) 

""nf [^J/ - €n - ^nf + t{Tnf + Iff + iD] = GnmTmf + Tf (^6^^ - t^^f^ - r„ {d^^^^ + t^^f) , 

where 

'Inn — '^^Gne^^Ge^^n + {GnkGkn/Pmk) , 

Cn = KP j de- GneG en/ {ei2 - e) + Im (GnkGkn/Pmk) ] 

7™j = TxhGfe^^Ge^^f + Re (GfkGkf/Pmk) , 

= KP j de- GfeGef/{eu -e)+lm (GfkGkf/Pmk) ; 

Ynf = T^^GneiiGeiif + R^ {GnkGkf/Pmk) , (3.36) 

5™^ = hV j de- GneGef/{ei2 - e) + Im (GnkGkf/Pmk) ; 

-fjl^ = TlhGfe^^Ge^^n + Re (GfkGkn/Pmk) , 

Sfn = hV J de- GfeG en/ iei2 - e) + Im (GfkGkn/Pmk) ] 
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Inn — 7r^Gne2^e2n + (GnkGkn/Pkn) , 

Kn = hV j de- GneGen/{e2 -e)+\m (GnkGkn/Pkn) ] 

-fjj = TihGfe^Gs^f + Re {GfkGkf/pkn) , 

^// = / ■ GfeGef/{e2 -e) + \m (GfkGkf/Pkn) ; 

'jUj = TlhGnciGeif + Re (GnkGkf/Pkn) , (3.37) 

^nf = KP j de- GneGef/{e2 -e) + \m (GnkGkf/Pkn) ; 

7/n = T^f^G ffr^Gir^n + Re {G fkGkn/Pkn) , 

S]n = KP j de- GfeGen/{e2 - £) + Im (GfkGkn/Pkn) ; 

7i{„ = TlhGnesGesn + Re {GnkGkn/Pkf) , 

= hV j de- GneGen/{e3 - s) + Im {GnkGkn/Pkf) ; 
^// = T^fiGfe.Ge.f + Re (GfkGkf/pkf) , 
Sjf = hV Jde- GfeGef/{es - e) + Im (G/fcCfc/Zp,./) ; 
-yl^f = nhGnssGesf + Re {GnkGkf/Pkf) , (3.38) 

= hV j de- GneGef/{e-i - e) + lm (GnkGkf/Pkf) ; 

7j„ = 7lhGfe.,Ge-,n + Re (GfkGkn/Pkf) , 

^fn = hV Jde- GfsGen/ies -e) + lm (GfkGkn/Pkf) - 



and 



£12 = h{uJi - UJ2), e2 = En - hu2, e3 = Ef - hws, (3.39) 
Pmfc = r^fc + z(u;i - u;2 - t^mfc), Pfcn = Tfcn + z(u;fc„ - cja), pfc/ = Tfc/ + z(cjfc/ - u;3). 

Besides the continuum states, a contribution of other non-resonant levels k is taken into 
account too and a sum over repeating k index is assumed. The contribution from these 
levels may occur commensurable to that of the continuum states. As seen from the equa- 
tions ()3.35p . the values 'jij and 6ij describe light-induced broadening and shifts of discrete 
resonances stipulated by the induced transitions between them through the continuum. 
The parameters qij = Sij/'jij are analogous to the Fano parameters for autoionizing states. 
Within the validity of (j3.3U|) . their dependence on the field intensities can be neglected. 
These parameters determine the most important features of the processes under inves- 
tigation since they characterize the relative integrated contribution of all off-resonant 
quantum states compared to the resonant ones. 
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With the aid of ()3.34|) . the equation for the dissociation rate fj3.32|) takes the form 



W = f, = 2 



7>n + 7^/^/ + 2Re (7/„r„/) . (3.40) 



The first two terms on RHS describe induced transitions to the continuum from the 
corresponding discrete levels, whose populations are determined by a variety of processes. 
Most important is the last term describing quantum interference, which occurs if two 
Lies overlap. Let us consider specific features of the processes in both open and closed 
energy level configurations. 

3.1.1 Open configuration 

In this case in the same way as above, one obtains 

rf{Tf + 27/;) = Qf- 2Re (r„,H7j„ + , (3.41) 

rniTn + 27;^„) = Qn - 2Re {iGnmrmn + r„/[(7/„ - i5/„]) • 
Then the solution of (j3.35j) and (j3.41|) can be presented in the form 

rnf = {(r„ - r„)|G„„r7r/(l - ^C/) - [rjlifil + ^qif) + r„,7::/l - tq:f)]Y]z"\ 
r„ = {L^C^ + L^S,)/{C,C2 - S.S^)^ r„ = {L^ + S^r^)/^, 

rf = [Qf + 2r^Bt + 2(r^ - rr,)Mt]Yf\ (3.42) 

The notations introduced here are as follows: 

Ci = Yr^^Yf - L,L2, C2 = [Tm + 2(0 - Mm]Yf - AM^Uf , 
Si = Km, + 2(G + Mr - D)]Yf + 2L2Mf, S2 = PYf + 2M/"Li, 
Fi = QmYf + 2QfMf, F2 = QnYf - QfL^, P = 2{G - + D), 
Li = 2(5^-^ - 2M^0, L2 = 2(5f - M/"), 17 = f ; - 2S|^, 
l^m = f „ + 2(6* - fir + -D),G= \Gran?Re {ymf/Y), 
5f =7L7LRe(6iV/Z), = |G_p7L7}/Re (6f /Z), 

D = \Gmn\Snnl7f^e [kr/{YZ)], J, 1 = 771,71, /; 

Y = ymnVrnf ' InnlTf^T"^ ^ ^ = Vrnfynfymn ' InnlTf^T"^ Vnf + \Gmn?y,nn. 
ymn ~ Tmn ~l" ^(^mn ^nn)' Z/m./ ~ ^mf ~l" ^(^m/ ^Jf)^ 

y^f = f „^ + t{n^j + 5L - h>l = h{l - tqjjil - iqlf), 
= Ul + V^J(1 + iqlf), bH = Kil - zqjjil + iq\^\ 

r^n = + 7^^, r„/ = + 7™ , r„ = r„ + 27I + Wn^, = + 27^^. (3.43) 
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First, consider the effect of the interference of two LICS in the case where the laser 
field at discrete transition is turned off {Gmn = 0). Then the solution ()3.42p takes the 
more simple form of 

Tn = {QnYf - 2Qnnn7j/Re (fe^Vz/n/) ) / A, 



rf = (QfYr, + 2g„7i7/"/Re (fof/y^^j )/A, 

r„ = f „ - 27L7;/Re (br/ynf) , 17 = f / - 27i7j/Re (fe^-Vl/n/^ 
A = YfY^ - 4(7l7i)2Re (^{Vz/n/) Re (b^^ /Vnf) ■ (3.44) 

The process of dissociation is described by the formula ()H.4()|1 with the aid of ()H.44|1 . All 
the terms proportional to Re (...) describe interference processes and disappear for either 
of the driving fields E2 or E3 being turned off. The corresponding spectral structures are 
power broadened (terms with 7) and of an asymmetric shape that is determined by the 
power-shifts (terms with 6) and by the Fano parameters q for the corresponding coupled 
excited levels. The position and shape of the laser-induced autoionizing-like resonance 
within the continuum can be varied with the change of the frequencies and intensities 
of the driving fields. Constructive interference can be turned into destructive by small 
variation of the overlap of two LICS, which forms a basis for coherent quantum control 
of photophysical processes by the driving lasers. This will be illustrated below with 
numerical examples. 

If either E2 or E3 is turned off (7;^^ = 0), the equations ()H.42|1 take the standard form 
that accounts for power-broadening of the level / or n depending on which field is turned 
off. When E^ = 0, we obtain 



\ /m 



2/L, 12 



mn I 



~l~ 2(rj7j -|- ) I C-mn I /Iz/mni 

QmJ^m/^n ~l~ "^i^Qm ~l~ Qn)|C'^mn| /Iz/? 



rj = Qf/T 



12 /u, 12 



, iy = 27lr„. (3.45) 
r +2fr +r -w )\G P/lu |2 

^m^'^y-'-m^^n "-'nm j \ mn\ / \ ymn\ 

These expressions describe a process of two-photon dissociation without any interference 
phenomena in the spectral continuum. Opportunities for manipulating dissociation and 
population transfer through the interference of two LICS follows from a comparison of 
equations (pH^ and 

3.1.2 Closed configuration 

Assuming f Wdt << 1 in (|3.29p . so that -|- r„ + ^ 1, the equations for levels 
populations in a closed scheme can be written as 

VnTn = -2Re [iGnmrmn + rnf{l]n + ^^/n)] + WnTm, 

r/f/ = -2Re r f^illf + i^if) + w jr^, = I - Tn - r f. (3.46) 
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Then, with the aid of ()3.35|) . the solution for the populations of the excited levels can be 
presented in the form 



Yf{wn + P) + Li{wf + 2M, 



A. 



nm 



Y„miwf + 2Mt^ 



Wf 



Ar\ 
- 2Mt") 



(3.47) 



(M;/-25f )(w7„ + P-Li). (3.48) 



The other notation are the same as in ()3.43|1 . Equations ()3.48|1 describe similar interference 
structures as discussed in the previous sections, but account for the effects of driving fields 
on incoherent excitation. For Ei turned off, these equations describing LICS reduce to 

Tn = {YfWn + KWf) /A'^, Tf = {Y^W f + kfWn) /A^, 

A'^ = {Yf + Wf){Yn + Wn)-knkf, (3.49) 

kn = 27l7j/Re (b^^/ynf) - kf = 27l7j/Re (^{Vyn/) - wj, (3.50) 

and the dissociation rate is calculated with the aid of equation (j3.4(Jj) . In the case of 
E3 = 0, dissociation from an incoherently populated level n is described by the equation 
in the common form of 



iy/2 = 7l^n,/(r„ + 27l + ^„). 



(3.51) 



Two-photon dissociation from the ground level, when the field E3 is off, can be analyzed 
with the formula 



W/2 



7;:„(w;n + 2|a 



Pr 

-L r 



i/\ymn\ ) 



(3.52) 



To a some extent, LICS can be treated as dressed multiphoton resonances accompanied 
by the interference of different quantum pathways. As follows from (|3.42j) and (|3.48p . when 
such resonances are not fulfilled and |2/n/|, |ym./| ^ l,\^rnn\, all interference structures 
disappear. The possibility of coherent quantum control of photodissociation based on 
the derived equations is illustrated below with the aid of a numerical model relevant to 
electronic-vibration-rotation terms of the molecule Na2 in the next section. 



3.2 Numerical analysis 

Many experiments on coherent control of branching chemical reactions have been carried 
out with sodium dimers Na2. In this case, level m of our model (Fig. ^ can be attributed 
to the state X^T,'^{v = 0, J = 45) for the close configuration and X^T,'^{v = 14, J = 45) 
for the open one, levels n and / to the states A^T,^{v = 6, J = 46) and B^Ilu{v = 
5, J = 45) correspondingly, and the continuum states e to the dissociation continuum 
Na(3s)+Na(3s). The relevant relaxation constants are: = 2 ■ 10'^ c~^, F„ = Fj = 
1.2 • 10® c~^. The relaxation rates for coherence (off-diagonal elements of the density 
matrix) are estimated as half of the sum of the corresponding rates for the diagonal 
elements. 
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3.2.1 Coherent control of populations and dissociation in folded schemes 

Numerical simulation for a closed-type scheme for the quasi-stationary case shows the 
possibilities of manipulation of continuum resonances through the effects of the interfer- 
ence caused by a strong laser field. Due to the action of strong radiation at a discrete 
transition, a split of level n appeares. Two LlCS in the dissociation continuum can be 
detected by variation of the detunings VLnf through the change of frequency of the laser 
field at the adjacent transition to the continuum {fe). 




Figure 10: Dissociation probability (a) and population of level / (b) as functions of 

^nf/^mf and Gmn/^mn- Here, 7nn/rm„, = 3, Iff/Tmf = 3, qnn = 0.2, Qff = -0.5, 
Qnf = 10, Wn/Tn = 0.1, Wf = 0, = 0. 

Figure ^1 shows that the dissociation rate can be considerably eliminated or alter- 
atively enhanced by varying either the intensity of the splitting field or the detuning Qnf- 
The spectral lineshape at constant field intensities has a well-known asymmetrical Fano 
profile. The plot (b) shows opportunities for manipulating the population of level /. 

The properties of LICS are essentially determined by the Fano parameters, which 
indicate an asymmetry in the distribution of the oscillator strength over the continuum, 
which is illustrated in Fig. El By varying the position of the LICS within the continuum, 
the intensities of the strong laser radiation and the one-photon detuning Qi, one can ma- 
nipulate the properties of the laser-induced structures. Large two-photon detunings lead 
to the vanishing of any interference phenomena in the continuum and thus is equivalent 
to the field being turned off. 
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Figure 11: Dissociation probability (a) and population of level / (b) as functions of 
^nf/^mf and of the Fano parameter Here, Gmnf^mn = 50. The other parameters 
are the same as in the previous figure. 




Figure 12: Dissociation probability (a) and population of level / (b) as functions of 
finf/^mf and of the Fano parameter Here, Gmn = 0. The other parameters are the 
same as in the previous figure. 

The above obtained results show that dissociation and population transfer through 
the continuum can be controlled by the field coupled to the discrete transition. Figure 
^1 displays a very trivial behavior for Ei turned off. Other opportunities of manipulating 
two-photon dissociation and populations transfer are associated with the control field 
coupling discrete and continuous states. 

In the next subsection, we investigate how such behavior may exhibit itself in difference- 
frequency four-wave mixing and nonlinear-optical generation of frequency-tunable radia- 
tion through the continuous states in folded schemes. 
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3.2.2 Coherent control of generated radiation in folded schemes 



Here we consider the case where the fundamental Ei and generated Es fields are weak. 
Then the solution of the density matrix equations f Section 13. 1|) for a closed-type scheme 
takes the form 



^mn mriUmf / iymrillmf ^1 ); 



(3.53) 



Similarly, for the case where the field Es is the probe radiation and Ei is the generated 
one, we obtain 



Rmn = iGmflfni^ - iq fn) / (ymnymf " 



(3.54) 



Then the wave equations for weak fields, coupled by the four-wave mixing process a;^ ^ 
Ui — UJ2 + ^3, have the form 

z) exp(iAA;z), 

dEi{z)/dz = i27ik[x?E2E;Esiz) exp{-iAkz), (3.55) 
where Ak = ks — ki + k2 — k^. Analogously to (|2.2p . the solution of equations ()3.55|) can 

(a) (b) 
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Figure 13: Dependence of the absorption coefficients ai/aio and as /aw, of the conversion 
rate b/ (aioa^o), and of the squared modulus of the nonlinear susceptibility Ix'-^-'P, on the 
detuning Qnf (a), and dependence of the quantum conversion efficiency r/q on the optical 
thickness and on Qnf (b). Here, C = 0.5, = 100, Qf = 6000, Qi/Tmn = 0, qjf = —0.5, 

Qnn = 0.2, Qnf = 0, (a,b). 

be presented like (I2.23|l with r/qg = Vqolx^^M'^ 9ngf, where gn = Inn/^mn, Qf = Jff/Tmf- 
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With the aid of the approximate expressions 



ttlO = 47lUJi\drnn\^/chTmn, "50 = 4:71UJ s\dmf\'^ / cHT 



IXol 



Tx/9>h) {1 + q^j)\dmndnedefdfm\ {'^mn'^ 



mf) 1 



(3.56) 



we obtain that -q^Q = (7r/2)^(l + g^^). Other notations are the same as in Section ITU 
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Figure 14: Dependence of the absorption coefficients ai/aio and as /aw, of the conversion 
rate fe/(aioa4o), and of the squared modulus of the nonhnear susceptibihty Ix'-^-'P on 
the detuning Qnf (a). Dependence of the quantum conversion efficiency ?7q on Qnf for 
zaio = 10^; zaio = 4 ■ 10^; zaio = 3 ■ 10^ (b). Dependence of the quantum conversion 
efficiency ?7q on the optical thickness and on Qnf (c). Here, gn = 6000, and the other 
parameters are the same as in the previous figure. 

Figure El illustrates the case where the strong fields are so large that power-broadening 
exceeds the half-width of the transitions {Tij), and also gf » gn- Then the generation 
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Figure 15: Dependencies of the absorption coefficients ai/aio and cts/aio, of the rate of 
conversion b, and of the square of the modulus of the nonhnear susceptibihty on 
the detuning Qnf (a). Dependence of the quantum conversion efficiency r/q on the optical 
thickness and on Qnf (b). Here, Qn = 10, Qf = 6.5. The other parameters are the same 
as in the previous figure. 

maximum falls to the parameter range for which absorbtion of the fundamental radiation 
becomes minimum and the rate of conversion is positive. With that, an oscillation regime 
appears. Within the interval of Qnf/^mf from —1000 to 3000, the rates of absorbtion and 
generation are almost equal {bo ~ 0), and the oscillations reduce to only one maximum 
of rjq. Figure El depicts the case where both driving parameters are equal {gn = Qf = 
6000). Here the maximum of rjq also falls within the range of minimum absorption of 
the fundamental radiation, and a conversion dip corresponds to the maximum of ai. 
Interference of LICS along with splitting of the discrete resonance give rise to complex 
structures in rjq while the waves propagate through the nonlinear medium [see Fig. EJc)]. 
For the driving parameters about three order less than in the previous case, and gn = gf, 
the implementation of interference allows one to achieve four orders of increase of the 
conversion rate bo (Fig. fT^ . The maximum in rjq becomes somehow smaller but is reached 
at a much shorter length of the nonlinear medium. With greater Fano parameters the 
potential enhancement of the output of generation through manipulating the interplay of 
LICS becomes greater. As shown in Fig. UHl ?7q — > 1 with g„/ — >• oo provided by the 
appropriate choice of the control fields intensity. Maxima of the quantum efficiency of 
conversion process reach 0.996 for g„j = 100, 0.919 for g„j = 5 and 0.569 for g„y = 0. 
At greater magnitudes of the resonance broadening increases, which leads to the 
requirement of longer nonlinear medium. This is seen from comparison of Figs. ITTT a) and 
ITTT b). It is seen from comparison of Figs. ITBT b) and llSf b) that most of the opportunities 
for manipulation appear at moderate detunings flmn- 
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Figure 16: Dependence of the quantum conversion efficiency rjq on the optical thickness 
and on Qnf for different Here, gn = 10, gf = 20; = 100 (a), g„/ = 5 (b), g„/ = 
(c). The other parameters are the same as in the previous figure. 
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Figure 17: Dependence of the absorption coefficients ai/aio and as /aw, of the rate of 
conversion 6, and of the squared modulus of the nonhnear susceptibihty on the 

detuning VLnj. Here, g„/ = 5 (a), g„j = (b). The other parameters are equal to those 
in Fig. [ini 



(a) (b) 




Figure 18: Dependence of the absorption indices ai/aio and as/aw, of the conversion 
rate 6, and of the squared modulus of the nonlinear susceptibility on the detuning 

Vtnf (a). Dependence of the quantum conversion efficiency ?7q on the optical thickness 
and on VLnf [max(?7g) = 0.87] (b). Here g„/ = 5, VLmn/'^mn = 100. The other parameters 
are the same as in the previous figure. 

Conclusions 

The technique of quantum control of such processes as ionization, chemical reactivity 
through dissociation of molecules and population transfer, and four-wave mixing is fur- 
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ther developed. Other applications involve discrete and broad-band spectra in solids. 
The approach is based on manipulating constructive and destructive interference of sev- 
eral quantum pathways involving both discrete energy levels (bound states) and energy 
continua (free states). Analytical and numerical solutions of the set of coupled density- 
matrix equations and Maxwell equations for travelling electromagnetic waves are found 
for the continuous-wave regime. 

Novel opportunities are shown for suppression or alternatively enhancement of pho- 
tophysical processes as well as photochemical processes related with branching chemical 
reactions. These opportunities are associated with the overlap of two laser-induced con- 
tinuum structures induced by two control fields, while the third strong field controls a 
discrete transition. The expressions obtained and numerical models developed are used to 
demonstrate the feasibility of manipulation of the spectral characteristics of the absorp- 
tion and dispersion both for discrete transitions and the spectral continua in the presence 
of additional strong laser radiations. Related opportunities to form laser-induced trans- 
parency and to enhance four-wave-mixing nonlinear-optical polarizations are investigated 
as well. It is shown that the implementation of such nonlinear interference effects makes it 
possible to utilize the advantages of multiple resonances and strong radiations for consider- 
able improvement of the generation of short-wavelength radiation. Among the important 
features discovered is the fact that the constructive or destructive nature of interference 
is governed not only by the ratio of the intensities and by the detunings from the discrete 
resonances, but also by the composite Fano parameters for transitions between high-lying 
discrete levels over the continuum states. The results can be generalized to higher-order 
processes. The matrix elements of the interaction Hamiltonian should then be replaced 
with the corresponding composite multiphoton matrix elements. 
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